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Abstract. We present observations of XRF 050406, the first burst detected by Swift showing a flare in its X-ray light curve. 
During this flare, which peaks at ~ 210 s after the BAT trigger, a flux variation of 6F/F ~ 6 in a very short time <5i/fpeak <K 1 
was observed. Its measured fluence in the 0.2-10 keV band was ~ 1.4 x 10~ 8 erg cnr 2 , which corresponds to 1-15% of the 
prompt fluence. We present indications of spectral variations during the flare. We argue that the producing mechanism is late 
internal shocks, which implies that the central engine is still active at 210 s, though with a reduced power with respect to the 
prompt emission. The X-ray light curve flattens to a very shallow slope with decay index of ~ 0.5 after ~ 4400 s, which also 
supports continued central engine activity at late times. This burst is classified as an X-ray flash, with a relatively low fluence 
(~ 10~ 7 erg cirr 2 in the 15-350 keV band, £ iso ~ 10 51 erg), a soft spectrum (photon index 2.65), no significant flux above ~ 50 
keV and a peak energy E p < 15 keV. XRF 050406 is one of the first examples of a well-studied X-ray light curve of an XRF. 
We show that the main afterglow characteristics are qualitatively similar to those of normal GRBs. In particular, X-ray flares 
superimposed on a power-law light curve have now been seen in both XRFs and GRBs. This indicates that a similar mechanism 
may be at work for both kinds of events. 
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1. Introduction 



The Swift Gamma-Ray Burst Explorer jGehrels etalj E004) 
was successfully launched on 2004 Nov 20. Its payload in- 
cludes one wide-field instrument, the gamma-ray (15-35 
keV) Burst Alert Telescope (BAT; Barthelm vetafl l2005al) . 
and two narrow-field instruments, the X-Ray Telescope (XRT; 
iBurrows et alj |20 05al) and t he Ultraviolet/Optical Telescope 
(UVOT: iRoming et aljEooll) . BAT detects the bursts, calcu- 
lates their position to ~ 1—4' accuracy and triggers an au- 
tonomous slew of the observatory to point the two narrow- 



field instruments. The XRT, which operates in the 0.2-10keV 
energy range, can provide ~ 5" positions, while the UVOT, 
which operates in the 1700-6000 A wavelength range, can fur- 
ther refine the afterglow localization to ~ Of.' 5. With its unique 
fast re-pointing capabilities Swift set out to investigate the very 
early phases of gamma-ray burst (GRB) afterglows, beginning 
as early as one minute after the BAT trigger. During the initial 
activation and calibration phases, which ended on 2005 Apr 5, 
Swift discovered 25 GRBs. The narrow-field instruments were 
re-pointed towards seven of them within a few hundred sec- 
onds, and such is the case for GRB 050406. 
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On 2005 Apr 6 at 15 :58:48.40 UT BAT triggered on GRB 
050406 (trigger 1 13872: |Parsons et alJl2005l) . and located it at 
RA(J2000) = 02 h 17 m 53 s , Dec(J2000) = -5£Ky52^witl^n 
uncertainty of 3 arcmin (95% containment; Kri mm et al.l2 005). 
The derived value for the time during which 90% of the burst 
fluence is observed was Tgo = 5+1 s in the 1 5-350 keV band. In 
the 15-25 keV band the light curve peak had a fast-rise, expo- 
nential decay (FRED) profile, while in the 25-50 keV band, the 
shap e was more sym metric, with the peak starting ~ 2 s earlier 
jRrimm et al.l f2005). Both the peak and time-averaged spectra 
were well fit by a simple power-law with a time-a veraged spec- 
trum photon index of 2.38+0.34 (90% confidence : lKrimm et alJ 
120051) . The fluence in the 15-350 keV band was 9.0 x 10" 8 erg 
cm' 2 . The gamma-ray characteristics of this burst, i.e. the soft- 
ness of the observed spectrum and the absence of significant 
emission above ~ 50 keV, classify GRB 050406 as an X-ray 
flash (XRF; iHeise et alJEool . From now on, we shall there- 
fore refer to this event as XRF 050406. 

Swift executed a prompt slew. The XRT imaged the BAT 
field only 84 s after the trigger but no bright X-ray source 
could be detected within the field of view. However, a refined 
on-gro und analysis revealed a previously uncatalog ued X-ray 
source (Cu sumano et al 1l2005aHCapalbi et al.l20 05). From the 
very first examination of the down-linked data it was clear that 
the afterglow of this burst was peculiar. Indeed, after an initial 
decay, the X-ray count rate bega n rising, peaking at « 220 s, 
and subsequently decaying again (Capal bi et al.l2005l) . 

Ground-based observations started as soon as the burst dis- 
covery was reported via the GCN network. The Magellan/Clay 
telescope imaged the XRT error circle with LDSS-3 in the 
R and i bands and found a single faint source (R = 22.0 ± 
0.09 mag, 7.8 hr after the burst) located at RA(J2000) = 
02 h 17 m 52 s :3, Dec(J2000) = -50° 1 1 ' 1 5" with an uncertainty 
of ~ (y.'5 in each coordinate jBerger et alJl2005albt) . Similarly 
to XRT, UVOT also imaged the field at the end of the slew 
(starting from ~ 88 s after the trigger) and thou gh it failed to 
detect the afterglow on-board ( Lands man et alJl2005t) . subse- 
quent on-ground analysis revealed a source within the XRT er- 
ror circle at the 4.3- (19.0 mag), 3.0- and 2.5-cr detection lev- 
els in the U, B and V bands, respectively. The UVOT position 
was RA(J2000) = 02 h 17 m 52 s :2, Dec(J2000) = -50°1 1'15'.'8, 
consistent with the Magellan one. By the time the second 
UVOT observation (1.3 hr later) was performed, the source 
was not detected in the U bancL confi rming it as the afterglow 
of XRF 050406. ISchadv et alJ ( l2006h obtained an estimate of 
Z = 2.44 ± 0.36 from fitting the broad band spectrum (com- 
bined UVOT and XRT data). 

In this paper we present observations of the first Swift burst 
where a flare is clearly detected in its X-ray light curve, during 
which the source count rate increased by a factor of ^ 6. This 
feature had never been observed before in Swift data, and had 
rarely been observed before in any X-ray afterglow JPiro et alJ 
2005). This paper is organized as follows. In Sect. |3 we de- 
scribe our observations and data reduction; in Sect. [3] we de- 
scribe our spatial, timing and spectral data analysis; in Sect. 
Uwe discuss our findings. Finally, in Sect. [5] we summarize 
our conclusions. Throughout this paper the quoted uncertain- 
ties are given at 90% confidence level for one interesting pa- 



rameter (i.e., Ay 2 = 2.71) unless otherwise stated. Times are 
referred to the BAT trigger To, t = T - To. The decay and 
spectral indices are parameterized as follows, F(v, f) oc t~ a v~P, 
where F y (erg crrT 2 s Hz -1 ) is the monochromatic flux as a 
function of time t and frequency v; we also use F = B + 1 as the 
photon index, N(E) oc E~ r (ph keV -1 crrT 2 s ). 

2. Observations and data reduction 

2.1. BAT observations 

Table ^reports the log of the observations that were used for 
this work. The BAT data were analyzed using the standard BAT 
analysis software distributed within FTOOLS v6.0. The burst is 
detected in the first two standard bands (15-25 and 25-50 keV) 
while virtually no flux is observed above 50 ke V. We find Tgo = 
6.1 + 1.0 s in the 15-150 keV band. 

The BAT spectra were extracted over the full time interval 
over which the burst was detected (r to t), in the interval covering 
the 1-s peak r pea k, and for the Tgo and Tso intervals. Response 
matrices were generated with the task batdrmgen using the 
latest spectral redistribution matrices. For our spectral fitting 
(XSPEC vl 1.3.2) we considered the 15-150 keV energy range. 
All spectra are well fit with a simple power law with photon in- 
dex T y ~ 2.65 (see details in Tabled). There is no evidence of a 
spectral break within the BAT energy range, thus constraining 
the peak energy E p < 15 keV. The indices are st eeper (softer) 
although consistent with the ones reported by iKrimm et alJ 
(2005), due to the different energy ranges used for the spec- 
tral fitting. No significant improvem ents are found usi ng either 
a cutoff power-law or a Band model IBand et alJl993l) . The 1-s 
peak photon flux was (2.3+ 28 )x 10~ 8 erg cm" 2 s" 1 (15-350keV 
band), while the fluence was T = (1.0*H?) x 10" 7 erg cm -2 
(15-350 keV band). This fluence corresponds to an isotropic- 
equivalent energy = (1.4+Jg) x 10 51 erg (in the rest fram e 
52-1204 keV) assuming z = 2.44 ± 0.36 JSchadv et all2006l) . 

2.2. XRT observations 

In order to cover the dynamic range and rapid variability ex- 
pected from GRB afterglows and to provide rapid-response, au- 
tomated observations, XRT was designed to support different 
readout modes that optimize the collected information as the 
flux of the burst diminishes. The XRT supports four major read- 
out modes, one imaging (IM), two timing, Piled-up/Low-rate 
Photodiode (PuPD and LrPD) and Windowed Timing (WT), 
and one Photon-Counting (P C). A deta iled description of XRT 
modes can be found in iHill et all J2004h . In the nominal oper- 
ating state the mode switching is based on the source flux and 
is fully automated (auto state) to minimize pile-up in the data. 

The XRT observations of XRF 050406 started on 2005 Apr 
6 at 16:00: 12 UT, only 84 s after the trigger, and ended on 2005 
Apr 22, thus summing up a total net exposure (in PC mode) of 
~ 163 ks spread over a ~16d baseline. The monitoring is or- 
ganized in 9 observations (000, 001, 002, 005, 006, 008, 009, 
010, 011) and 183 snapshots (continuous pointings at the tar- 
get). This was the first burst to occur after the formal end of 
the calibration phase (2005 Apr 5), and the first (000) obser- 
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Fig. 1. XRT image of XRF 050406, obtained from the to- 
tal ~163ks PC mode data. The field is centred on the 3' 
radius BAT error circle. Also sho wn is the XRT 4'.'2 er- 
ror circ le, as well as th e Magellan (iB erger et alJl2005al) and 
UVOT dRol et alj E005) optical counterpart positions; the op- 
tical points are so close they cannot be distinguished on this 
scale. S2 is a serendipitous source located at RA(J2000) = 
02 h 17 m 52. s 9, Dec(J2000) = -50°10'36'.'l. 

vation was performed as an automated target (AT) with XRT 
in auto state. Therefore, during observation 000 the automated 
mode switching made XRT take an initial 2.5 s image (IM at 
t = 84 s), immediately followed by one PuPD (f = 90 s) and 
one LrPD (? = 91 s) frame. Then at t — 92 s a series of 5 WT 
frames was taken until the on-board measured count rate was 
low enough for XRT to switch to PC mode (t = 99 s). After 
this, XRT repeatedly switched between WT and PC modes be- 
cause of an increased background level (see below). Since the 
signal-to-noise (S/N) in these late WT frames is low, we did 
not include them in our analysis (TableQ. 

The XRT data were first processed by the Swift Data Center 
at NASA/GSFC into Level 1 products (calibrated and quality- 
flagged event lists). Then they were further processed with the 
XRTDAS (v 1.4.0) software package written by the Agenzia 
Spaziale Italiana (ASI) Science Data Center and distributed 
within FTOOLS v6.0 to produce the final cleaned event lists. 
We ran the task xrtpipeline (v0.8.8) applying standard fil- 
tering and screening criteria, i.e., we cut out temporal intervals 
during which the CCD temperature was higher than -47 °C, 
and we removed hot and flickering pixels. These are present 
because the CCD is operating at a temperature higher than the 
design temperature of -100 °C due to a failure in the active 
cooling system. An on-board event threshold of ~0.2 keV (un- 
reconstructed pulse-height PHAS[1]> 80) was also applied to 
the central pixel, which has been proven to reduce most of the 
background due to either the bright Earth limb or the CCD dark 
current (which depends on the CCD temperature). These two 



sources of background are the main reason for the switching 
between PC and WT mode even when the source count rate is 
below 1 counts s _1 . 

Throughout the monitoring campaign the CCD tempera- 
ture was < -50 °C, with the exception of part of observations 
002 and 005, where it became as high as -43.5 and -45 °C, 
respectively; those data were therefore screened out. For our 
analysis we further selected XRT grades 0-12 and 0-2 for PC 
and WT data, respe ctively (according to Swift nomenclature; 
burrows et all2005al) . 

3. Data analysis 

3.1. Spatial analysis 

Figure^shows the 163 ks XRT image accumulated in PC mode 
in the 0.2-10 keV energy band. We detected two previously 
uncatalogued sources within 1 arcmin of the optical burst co- 
ordinates. The brightest uncatalogued source, which we iden- 
tified as the fading X-ray counterpart of the burst, is present 
in the first four XRT snapshots. The source is piled-up during 
the initial 500 s of PC data. Therefore, to obtain an unbiased 
position, we rely on the remainder of the PC data in the first 
observation, which has a net exposure of 49.8 ks. We used the 
xrtcentroid task (vO.2.7) and found that the afterglow posi- 
tion is RA(J2000) = 2 h 17 m 52 s .4, Dec(J2000) = -50° 1 1'13'.'6. 
We estimate its uncertainty to be 4'.'2 (90% confidence level). 
This position takes into account the correction for the mis- 
alig nment between the telescope and the satellite optical axis 
(Mor ettiet J1 I2005). Figure ^ s hows the XRT erro r circle, as 
well as the 3' BAT error circle (Krim m et alJl20 05: 95% con- 
tainment) and the optical c ounterpart coordinates determined 
by M agellan dBereeret al]l2005al) and by UVOT dRol et alJ 
2005). The XRT coordinates are 23" from the BAT ones, and 
l'/6 and 2'.'8 from the Magellan and UVOT ones, respectively. 
XRF 050406 was detected (XIMAGE v4.3) in the first four snap- 
shots individually, but not from the 5th on. The second source, 
S2, is located at RA(J2000) = 02 h 17 m 52 s :9, Dec(J2000) = 
-50°10'36'.'l and has a constant rate (3.8 + 0.7) x 10~ 4 counts 
s _1 throughout the observation campaign. 

3.2. Temporal analysis 

During the first 500 s of the XRT observation the intensity of 
the afterglow was high enough to cause pile-up in the PC mode 
data. To account for this effect we extracted the source events 
within an annulus with a 30-pixel outer radius (~ 71") and a 2- 
pixel inner radius. These values were derived by comparing the 
observed and nominal PSF For the PC data collected after the 
first 500 s, the entire circular region (30-pixel radius) was used, 
instead. In both cases we further disregarded data within a cir- 
cular region centred on the serendipitous source S2 (which lies 
within the 30-pixel PC source extraction region) with a 7.17 
pixel (17") radius. The WT data were extracted in a rectangu- 
lar region 3 1 pixels long along the image strip (and 20 pixels 
wide), which excludes the data from the source S2. The se- 
lected extraction regions correspond to ~ 69 % (piled-up PC), 
~ 95 % (non piled-up PC), and ~ 94 % (WT) of the XRT PSF. 



P. Romano: X-ray flare in XRF 050406 



CO 
CO 

C 

o 

o 

CD 

co 
Si 

c 

3 

O 

o 



10.0000 r 



1.0000 r 



0.1000 r 



0.0100 r 



0.0010 r 



0.0001 



~ i — i — i — i 1 1 1 1 



1 1 1 1 1 1 — i — i — i 1 1 1 1 1 




1000 



_l I I 



_l I I 



_l I I 



_l I I 



10' 



10° 10* 10° 

Time since BAT trigger (s) 



10 



6 



Fig. 2. X-ray light curve of the XRF 050406 afterglow in the 0.2-10ke V energy band. The curve is background-subtracted and 
the time is referred to the BAT trigger, 2005 Apr 06 at 15:58:48.4 UT JParsons et alJl2005l) . The last point after 10 6 s is a 3-cr 
upper limit. Inset: details of the first ~ 1000 s, which include data in all XRT modes. The (yellow) diamonds represent LrPD 
mode data taken during the latter portion of the slewing phase; the (cyan) triangle is the initial IM point (84 s after the trigger, 
see Table the downward-pointing arrow is a LrPD limit (pointing, 91 s after the trigger), the (blue) circles are WT mode 
data (starting from 92 s after the trigger), and the (red) squares are PC mode data (starting from 99 s after the trigger). The data 
have been corrected for pile-up (where appropriate) and PSF losses. The solid (red) line represents the best-fit broken power-law 
model to the light curve (excluding the flare). 



To account for the background, data were also extracted in PC 
mode within a circular region (radius 130" = 54.8 pixels) and 
in WT mode within a rectangular box (40x20 pixels), in loca- 
tions far from background sources. The mean PC background 
in the 0.2-10 keV band was found to be constant throughout 
the observations and, normalized to the PC source extraction 
region, it had a value of ~ 2.6 x 10~ 3 counts s _1 . The mean WT 
background in the same energy band and normalized to the WT 
source region was ~ 4.6 x 10 -2 counts s . 

Figure |2 shows the background-subtracted light curve ex- 
tracted in the 0.2-10 keV energy band, with the BAT trigger 
as origin of time. We considered WT data for the first snap- 
shot of the first observation, and PC data for all 9 available 
observations (see Table 0. During the initial phases of the af- 
terglow evolution it < 4 x 10 4 s) we binned the source counts 
with a minimum of 30 counts per time bin, and dynamically 
subtracted the normalized background counts in each bin. The 



PC mode data were corrected for the effects of pile-up. We 
note that, by keeping to the minimum number of counts per 
time bin criterion, we created several bins during the first snap- 
shot, but subsequently needed to merge data belonging to snap- 
shots 1 and 2 (point at ~ 4ks), then from snapshots 3 and 4 
(point at ~ 20 ks), and later on from snapshots 5 through 8 
(point at ~ 35 ks). Afterwards, we used XIMAGE with the op- 
tion S0STA, which calculates vignetting- and PSF-corrected 
count rates within a specified box, and the background in a 
user-specified region. To ensure uniformity with the early light 
curve, the background was estimated in the same region as the 
one used for the initial part of the light curve. We thus obtained 
a signal-to-noise ratio S/N <: 3 (the only exception being the 
point at ~ 33 ks which has S/N <: 2). The last point is a 3-cr 
upper limit. This latter method is preferred for the construction 
of the late part of the light curve since it better accounts for 
the background in a low-counts regime. We note, however, that 
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extracting the light curve in the same 30-pixel source region 
up to the end of the last observation, we obtained fully con- 
sistent results, albeit with a noisier light curve. We also note 
that the residual contribution of the serendipitous source S2 
within the source extraction region is 5= 19% of the S2 counts, 
which corresponds to ^ (7 + 1) x 10~ 5 counts s" 1 . Therefore, 
S2 only makes a marginal contribution to the afterglow light 
curve, which amounts to <20% of the last point. 

The light curve clearly shows a complex behaviour, with 
a power law decay underlying a remarkable flare which peaks 
at » 210 s after the BAT trigger (see Fig. inset). To fit the 
light curve we used the BAT trigger as reference time and we 
only considered spectroscopic-mode data obtained while XRT 
was pointing, thus excluding the early LrPD, the LrPD upper 
limit and the IM point. Further excluding the data taken during 
the flare (180s < t < 300 s), a fit with a simple power law 
yields - 4.32 (12 degrees of freedom, d.o.f.), which is 
unacceptable. A fit with a broken power law F (f) = Kr" 1 for 
t < f b and F(t) = Kt b ai (t/t b y a2 for t > f b , where t b is the time 
of the break, yields ct\ = 1 .58+g }| and a 2 = 0.50 ± 0.14, and a 
break at ~ 4200 s after the BAT trigger. This latter model yields 
a good fit (xLj, = 1.20, 10 d.o.f.), a significant improvement 
over the simple power law (null hypothesis probability = 1 .7 x 
10" 3 , equivalent to 3.2 &), but some of the parameters are not 
well constrained. Alternatively, a fit with two smoothly joined 
power laws F(f) = K'[(t/t b y a < + (f/? b )" ff2 ] yields X 2 kA = 1-29 
(10 d.o.f.) with similar values for the inferred parameters. A 
summary of the fits to the light curve can be found in Table 
As a reference, the 0.2-10 keV unabsorbed flux at fb is (4 ± 
1) x 10" 13 erg cm" 2 s" 1 (we adopted a count rate to unabsorbed 
flux conversion factor of 6.5 x 10" 11 erg cm" 2 count" 1 , obtained 
from the best fit models derived in Sect. 13. 3> and the luminosity 
in the 0.7-34.4 keV band is (1.9 ± 0.9) x 10 46 erg s _1 . 

During the flare a rebrightening of the source by a factor 
of ^ 6 in flux was observed between t ~ 154 s and the peak 
at ~ 210 s. Both the rising and the falling part of the flare had 
very steep slopes that, when fit with a simple power law, yield 
£*i,flare = — 5.8*i"i an d a 2,flare = 6.7 + 1.0. When the underly- 
ing power-law afterglow is subtracted, the fit yields ari >nare = 
—6.8*2i m ^ a 2,flare = 6.8*|q and the peak is at 213 + 7 s from 
the BAT trigger. In all cases the errors are dominated by the un- 
certainty in the placement of the flare boundaries. The flare can 
also be characterised, as a simple parametric description, as a 
Gaussian line. A combined broken power law and Gaussian 
model fit yields a peak at 21 1.1 s (61.4+JoS in the rest 
frame) and a width 17.9±' 2 6 3 s 0r 2 d = 1.58, 17 d.o.f.). In this 
case the ratio of the characteristic time-scale and the peak time 
is <Sf/fpeak ~ 0.08 or 0.20, when using the Gaussian width or its 
FWHM (42.2+ 2 q " s), respectively. In either case, (Jf/f pea k « 1, 
which puts severe constraints on the emission mechanisms that 
can produce the flare. We shall address this issue in the discus- 
sion section. Integration of the Gaussian best-fitting function 
yields an estimate of the fluence of the flare, (1.4 + 1.0) x 10" 8 
erg cm" 2 , corresponding to an energy of (2.0 ± 1.4) x 10 50 erg. 
The large error reflects the uncertainty on the actual model used 
for the integration of the flare. 

We also extracted events from the first snapshot WT data 
in two more energy bands, 0.2-1 keV (soft, S) and 1-10 keV 



1 

(T) 0.2-10 keV ^ j 


(a) : 






•nm|iim|mi 

I (S) 0.2-1 keV _j_ 

\ 


(b) i 




,,,,,,,,, 
L (H) 1-10 keV 

\ — *— -t- " + "-+-^ 

\ ' 1 ' ' 1 ' 1 ' ' 1 ' 1 1 ' 


< < < 1 < 

(c) 1_ 

1 , , , , 1 , ,; 


L. .+ .+!-+: 


1 . . 


(d) j 

, 1 , , : 



100 200 300 400 

Time since BAT trigger (s) 



Fig. 3. WT background-subtracted light curves, (a): Total band 
(T, 0.2-10 keV). (b): Soft band (S, 0.2-1 keV). (c): Hard band 
(H, 1-10 keV). (d): Ratio of hard to soft count rates. 

(hard, H), as well as the total band, 0.2-10 keV. We used the 
same regions as the ones described above, a constant time bin- 
ning of 30 s and dynamically subtracted their respective back- 
grounds. Figure[3]shows the three background-subtracted light 
curves, as well as the ratio H/S. Indeed, during the rising por- 
tion of the flare the hard band flux increases by a factor of ^ 6 
while the soft band flux only increases slightly, so that the spec- 
trum of the flare starts off harder than the underlying afterglow, 
and then evolves into a softer state as its flux decreases; this 
can be seen in the following time bin, when the soft band flux 
peaks with a flare to pre-flare flux ratio of ~ 3.5. This yields 
an indication of spectral evolution during the flare as a ~ 3-<x 
excess over a constant fit to H/S. It should be noted that this 
beha viour is reminisc ent of that observed in the prompt emis- 
sion ( For d et al ■119951) . with the harder band peak preceding the 
softer band peak. 

At t ~ 1.7 x 10 5 s a second faint bump is observed. Its 
significance is not high, since it is detected as a 2-cr excess over 
the underlying afterglow. Similar late-time bumps have been 
observed in other S wift-detected GRBs (e.g. GRB 050502B; 
iFalconeet all2006l) . 

3.3. Spectral analysis 

The afterglow of XRF 050406 was very faint, hence it is not 
possible to perform time-resolved spectroscopy to distinguish 
the spectral properties of the afterglow proper from the ones 
of the flare observed in the light curve. Therefore, we pro- 
ceeded as follows. Spectra of the source and background were 
extracted in the regions described in 33. 21 from the first obser- 
vation (000) event files. PC and WT spectra were extracted dur- 
ing the first ~ 500 s of the PC observation (see Table|2]for times 
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referred to To), when PC data are piled-up and when the flare 
is observed in the light curve. We also extracted PC spectra 
after the first 500 s during the first 4 snapshots. For the latter 
we used a circular region with a 10 pixel radius (correspond- 
ing to ~ 80 % of the XRT PSF) to mi nimize the b ackground 
and to be able to use the Cash statistics JCashl 19791) . Ancillary 
response files (ARF) were generated with the task xrtmkarf 
within FTOOLS v6.0 using the latest ARF distribution (v003). 
These ARFs account for different extraction regions and PSF 
corrections. We used the latest spectral redistribution matrices 
(v007). The energy ranges adopted for spectral fitting were 0.5- 
10 keV and 0.2-10 keV for WT and PC, respectively. 

We first performed a fit with an absorbed (wabs in XSPEC) 
power law to the WT data (166 counts), which were rebinned 
with a minimum of 20 counts per energy bin to allow x 2 fit- 
ting within XSPEC. The Hydrogen column was initially kept 
as a free parameter, an d then frozen to the Gala ctic value 
(N% = 2.8 x 10 20 cm' 2 , bickev & Lockman1ll990l) when the 
fit yielded a value lower than (although consistent with) the 
Galactic one. The fit was good, x 2 ^ = 10 for 6 d.o.f., and 
yielded F = 2.11*!H!i. We then performed a fit with the same 
model to the remainder of the PC data during snapshots 1 
through 4 using Cash statistics which is more appropriate given 
the low number of counts (21 un-binned counts) and calculated 
the goodness of the fit via 10 4 Montecarlo simulations. The fit 
was good and yielded consistent results. We also performed si- 
multaneous fits to the WT and PC (60 counts) spectra extracted 
during the first ~ 500 s (using x 1 statistics) and of the PC data 
alone (using Cash statistics), also obtaining consistent results. 
Table |2 summarizes the results of the fits. We note that, given 
the current goodness of the XRT calibration (5% systematic 
uncertainty for all observing modes and grade selections in the 
0.5-10 keV range; e.g., Rom ano et aDf2 005). an excess of Nu 
cannot be excluded and we find a 3-cr upper limit to the to- 
tal (Galactic plus intrinsic) Hydrogen column along the line of 
sight of jVh < 9 x 10 20 cirT 2 . 

We can therefore conclude that, during the first 600 s after 
the burst, which include the X-ray flare observed in the light 
curve, the mean photon index is T = 2.1 +0.3, and that the pho- 
ton index does not vary after the end of the flare. However, we 
do have clues regarding the presence of spectral evolution dur- 
ing the flare coming from the hardness ratio analysis (Sect. 13. 2\ . 
even though the statistics are not high enough to show it in 
the spectral analysis. As we will discuss later (Sect. l4~2"l . other 
afterglows with larger amplitude X-ray flares demonstrate a 
strong spectral evolution of the flares. 

4. Discussion 

4.1. Gamma-ray properties: similarity of XRFs and 
GRBs 

The duration of this burst (Tgo = 5 + 1 s in the 15-350 keV 
band) pla ces this burst in the sho rt tail of the long GRB pop- 
ulation (IKouveliotou et al.lll993l) . Its fluence is relatively low 
(1.0 x 10" 7 erg cm -2 in the 15-350 keV band) but not un- 
usually faint. The gamma-ray characteristics of t his burst are 
consistent with a classification as an X-ray flash faeise et alJ 



2001!), or as an "X-ray rich GRB" (XRR). The softness of the 
observed spectrum, which is well fit in the 15-150 keV band 
with a simple power law with photon index F y = 2.65, and with 
no significant emission observed above ~ 50 keV, implies that 
the peak energy is below the BAT bandpass (£„ < 15 keV). T he 
operational definition of XRFs/XRRs (e.g. lLamb et al . 2004) is 
of a fast transient X-ray source characterized by a softness ra- 
tio SR= log[:F(2-30keV)/:F(30-400keV)] > for an XRF 
and -0.5 <SR< for an XRR. Extrapolation of the BAT spec- 
trum, with the assumption of E p < 2 keV, yields SR= 0.8*K, 
which classifies this burst as an XRF. However, a break in the 
spectrum may well be present in the 2-15 keV band. In the 
most conservative case, i.e. assuming no flux below 15 keV, 
this event would be an XRR GRB, with SR= -0.2+° \. 

The isotropic-equivalent gamma-ray energy of this event 
is £ iso = (1.4+^) x 10 51 erg (Sect. |2j), and this effec- 
tively _j2uts_XRF 050406 in the low-energy tail of GRB ener- 
gies ( Bloom et al.l|2003bl) . Assuming that the Amati relation 
llAmati et al J 120021) holds, we can infer a rest- frame E l p A ~ 
55 keV, which corresponds to an observer- frame E p ~ 15 keV. 
This value is consistent with the nondetection of E p in the BAT 
energy range. 

To date, X-ray afterglows of XRFs have b een detected 
in just a few cases (XRF 01 1030, XRF 20427: iBloom et all 
20033 iLevan et alJ 120051 XRF 03072 3: butler et alJ 120041 
XRF 040701: lFoxll2004t XRF 050315: Ivaughan et al.l l2006). 
This is one of the first examples of a well-studied X-ray light 
curve of an XRF. Its main characteristi cs are not qualitatively 
different from those of normal GRBs (iChincarini et al.l l2005: 
iNousek et alJl2006l) . As observations accumulate, it is becom- 
ing clear that these two classes of phenomena share many 
propertie s, and both have a fterglows with similar characteris- 
tics ( Sakamot o et al.lEo05l) . This is a clue that both types of 
events may have a common origin and is supported by re- 
cent evidence that some XRFs are associated with supernova e 
dSoderberg et al 120041: ISersier et al.l2005llFvnbo et alJ2004 . 



4.2. X-ray flares: evidence for prolonged engine 
activity 

The general behaviour of the afterglow of XRF 050406 is a 
typical one. The observed X-ray pho ton index (Fx = 2.1) 
is common among X-r ay afterglows ( Chincar ini et alJ 12005: 
be Pasauale et alJl2005l) . The light curve shows a break from 
a relatively steep decay (a\ = 1.58) to a flatter one (ai - 0.50). 
Its overall shape is simi lar to the one typically ob served by the 
XRT dChincarini et alJl2005l Rousek et alJl2006l) . even though 
the initial slope is less steep than average. 

The most striking characteristic of this burst is the strong 
flare in its X-ray light curve, a feature which had never been 
detected by Swift before and h ad been previous ly observed in 
very few GRBs ( GRB 970508.IPiro et alJll999l: GRB 011121 
and GRB 01 121 l. lPiro et alJl2005l) . The fluence of the flare is 
~ 1 .4 x 10~ 8 erg cm' 2 in the 0.2-10 keV band, which amounts 
to ~ 14 % of the observed (15-350 keV band) prompt fluence. 
A better estimate of the flare-to-prompt energy ratio would re- 
quire the knowledge of the prompt spectral energy distribu- 
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tion (SED). Since the actual peak energy of the prompt SED 
is unknown (E p < 15keV), the extrapolation of the BAT flu- 
ence to the XRT band is highly uncertain. For plausible val- 
ues of £p, the flare to prompt fiuence ratio is in the 1-10% 
range. The observed rebrightening is by a factor of 6 in flux, 
presents a peak at f pea k = 213 + 7 s and takes place on a very 
short timescale, with a ratio of the characteristic time-scale and 
the peak time <5f/f pea k ^ E Both the rising and the falling 
parts of the afterglow-subtracted flare had very steep slopes, 
cci.flare ~ —7 and a2,flare ~ 7, assuming the burst trigger as the 
time origin. 

According to the standard relativistic fireball model, the 
prompt emission is caused by internal shocks within the ex- 
panding fireball, while the afterglow is produced b y the fireball 
shocking the external me dium (external shocks. Ipiranl Fl 999 : 
IZhang & Meszarosl l2004l) . Available models to exp lain flares 
include refreshed shocks jRees and Meszaros| |l998|). external 
shocks with a clumpy medium (Lazzati et al. 200 21) and an- 
gular inhomogene ities in the outflow jFenimore et alj[l999: 
Nakarej: al]|2QQ3| ). However, it can be argue d jBurrowsetall 
2005b UZhang et alJ2005llNousek et alJ2006l) that such models 
cannot produce the observed large flux vari ations 6F/F„^ » 
1 in such short timescales 5f/f pea k ^ 1 dloka etalJ 12005). 
Similarly, none of the above mechanisms would explain the 
steep slopes observed in the flare. External reverse shocks, cre- 
ated when the fireball slows down because of the interaction 
with the external medium, are expected to emerge at optical 
and radio wavelengths, hence synchrotron self-Compton (SSC) 
must be invoked to produce emission in the X-ray band. This 
would require carefully balanced conditions ( Koba vashi et alJ 
l2005h. 

IPiroet al1 J2005) suggested that the X-ray flares observed 
in GRB 01 1 121 and GRB 011211 were due to the onset of the 
afterglow. The steep slopes and the short timescale variabil- 
ity can only be accounted for within th e thick shell scenario 
JSari & PirarJl999l) . lGalli & Pirol J2005I) successfully modeled 
XRF 011030 using this model. In this scenario, the emission 
before and after the flare is due to different processes (prompt 
tail and afterglow, respectively), hence a discontinuity in the 
light curve is generally expected underlying the flare. This is 
not the case for XRF 050406, where the same component de- 
scribes the X-ray emission both before and after the flare. Even 
if a fine-tuning may explain this particular event, the lack of 
a light curve break is comm on to a large fraction of the flares 
observed by Swift feurrows et alJ2005d) . Therefore, while the 
explanation of flares in terms of the afterglow onset is attrac- 
tive, it is unlikely to be applicable to the vast majority of the 
X-ray flares seen by XRT. 

A p romising mechanism to produce the flare is late internal 
shocks torn & Weill2005l IZhang etaljEool: iKing et aljEool 
IPerna et alJl2006h . which implies that the central engine is still 
active at t - 213 s, even though the prompt emission ended 
after t ~ 6 s. The late-time activity in this case must have 
a reduced power with respect to the prompt emission, as the 
relative fluences indicate. Such a mechanism would naturally 
explain the steep rise and decay slopes. Second, the energy re- 
quired to power the flare wou ld be much lower than in the other 
scenarios (Zhan g et al.l l2005). The indications of spectral evo- 



lution throughout the flare further support this interpretation. 
The flare appears to be harder than the underlying afterglow, 
which suggests a distinct origin for this emission. Furthermore, 
there are indications of spectral evolution, which shows the typ- 
ical hard-to-soft pattern. Such a behaviour is commonly ob- 
served in the prompt emission spikes of GRBs (e.g. lFord et alJ 
1995), which are produced in internal shocks. Further evidence 
of late engine activity comes from both the flat part of the light 
curve («2 ~ 0.5, see Sect. 14. 31 and possibly by the presence of 
the late-time bump observed at t ~ 1.7 x 10 5 s. 

Following the discovery of a flare in the af t erglow of XRF 
050406, initially reported by iBurrows et alJ J2005bl) , man y 
others were ide ntified: GRB 050502B JFalcone eTall 12006). 
GRB 050724 JBarthelmv et alJ l2005bh and GRB 050904 
JCusumano et alJ l2005<J) . just to mention a few. At the time 
of writing (2005 Oct), ~ 50 % of the bursts detected by XRT 
which were immediately re-pointed towards showed flares, 
making flaring quite a common behaviour. Furthermore, all the 
characteristics of the XRF 05040 6 flares have now bee n ob- 
served in most flaring GRBs (see Bu rrows et al.ll20 05c for a 
recent review). For example, highly significant spectral evolu- 
tion throughout the flare has been reported in GRB 050502B 
(which was the brightest observed so far) and GRB 050724. In 
several cases the flares present large amplitudes and occur on 
short timescales. Furthermore, several flares are often observed 
in the same event, at times ranging from ~ 100 s to 1 4 — 1 5 s 
after the burst. Finally, in most cases the afterglow is clearly 
present before the onset of the flare, and has consistent decay 
slope and flux levels with after the flare. The present case shows 
that flares are present both in XRFs and in GRBs. Since flares 
are likely tied to the central engine activity, this finding further 
supports the id ea that a similar m echanism is at work for both 
kind of events jFan & Weil2005l) . 

4.3. The X-ray afterglow light curve 

The prompt reaction of Swift has allowed us to observe the X- 
ray light curves of GRB afterglows starting from a few tens 
of seconds after the burst explosion. In most cases the X-ray 
light curves are characterized by an initial steep decay (up to 
~ 500 s) followed by a shallow decay, and then by a steeper 
decay with a seco nd break normally occurring at a few thou - 
sand seconds later JChincarini et alJ2005UNousek et alJ20 06). 
The early steep decay seen in the X-ray light curve can be 
explaine d as the tail of the prompt emission (however, see 
IPanaitescu et alJl2005h . The few cases where the XRT light 
curve lies well above the extrapolation of the prompt emission 
into the X-ray band can be explained either by a strong spec- 
tral evolution or by an X-ra y flare with the maxim um located 
before the XRT observation (Tagliafe rri et all200 5). There are 
othe r instances where the fi rst steep decay is not observed at all 
fe.g. lCampana et all2005l) . 

In the case of XRF 050406, however, the initial slope is 
shallower than the stee p values 3 ^ a ^ 5 observed in other 
early afterglows ( Taglia ferri et aT]l2005l). Moreover, the curva- 
ture r elation a = 8 + 2 JKumar & Panaitescu! l2000t iDermeil 
2004) is not satisfied, even after taking into account the ef- 
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fects pointed out bv lZhang et ail d2005t) that would alter such 
relation. Therefore, we also investigate whether the initial de- 
cline seen in XRF 050406 is consistent with afterglow emis- 
sion. Comparison of spectral indices and temporal decay slopes 
with theoretical r elativistic fireball models (e.g. Table 2 in 
IZhang etalJl2005h indicates that the first decay index ct\ - 
1.58 +0.17 and energy index/? = 1.1 +0.3 rule out fast cooling 
models (for which the injection frequency v m exceeds the cool- 
ing frequency v c ) for v < v m . For v > v m , the a(J3) = (3/3 -1)12 
closure relation is satisfied within the errors and an electron 
power-law distribution index p a; 2.5 is obtained. The same re- 
lation holds for the slow cooling regime (where v c > v m ) for 
v > v c (both wind and ISM). In this case a consistent solution 
is also found for v m < v < v c , although with a large p « 3. The 
ISM environment is favoured on the basis of a better satisfied 
closure relation. In conclusion, the spectral indices and tem- 
poral decay slopes of the first part of the X-ray curve can be 
interpreted in terms of relativistic fireball models, even though 
the large uncertainties associated with the slopes do not allow 
us to choose among the available models. 

An alternative explanation for the initial XRT emission is 
the presence of an additional flare which started before the be- 
ginning of the XRT observation, and of which we only see the 
decaying part. The superposition of two (and possibly more, 
fainter) flares would then mimic the initial steep power law de- 
cay. However, this interpretation seems less likely since recent 
Swift observations of X-ray flares within the first several hun- 
dred seconds of the prompt emission all had temporal decay 
indices much steeper than the observed XRF 050406 pre-flare 
index. 

At t ~ 4400 s the XRT light curve breaks to a 2 ~ 0.5. 
Such a flat decay cannot be explained in terms of the standard 
afterglow model. The only possibility would be to observe, in 
the fast cooling regime, the segment with v c < v < v m (where 
a = 0.25 is expected, marginally consistent with the observed 
value). However, the fast cooling regime is expected to end 
much earlier. To maintain the observed decay unbroken up to 
~ 10 6 s, large values of the equipartition parameters e e and sb 
or of the Compton parameter would be required. We consider 
this possibility quite unlikely. Another possibility is that the 
angular energy profile of the fireball is not trivial (a structured 
jet), so that emission coming from the (brighter) wings of the 
jet may in crease the observed flux as the fireball Lorentz factor 
decreases dPanaitescu et all2005l) . 

An interesting explanation for the shallow-decay phase is 
injection of new energy into the fireball through refre shed 
shocks dSari & MeszaroslkOOOt Ehang & Meszarosl200ll) . For 
this to happen, the energy release inside refreshed shocks must 
be sizeable, since the whole fireball dynamics has to be mod- 
ified. Assuming an energy injection rate E oc r q , we find q in 
the ra nge to 0.5 depending on the model details (Zhang ~t alJ 
2005). In this model, the initial part of the XRT afterglow light 
curve can be due to standard afterglow emission only if the fire- 
ball evolution is not influenced at these stages. Indeed, the en- 
ergy supply provided by refreshed shocks is steadily growing, 
and at the beginning it cannot alter the fireball dynamics. In this 
case, the break would identify the time when the new, injected 
energy is comparable to the fireball energy. On the contrary, if 



the first XRT phase were due to late engine activity, then the en- 
ergy injection could have begun much earlier and its emission 
would have been masked. 

Integration of the light curve from the onset of the flat slope 
phase yields T * 3 x 10~ 8 (f en d/7.6 x 10 5 s) 5 erg cirT 2 , where 
/ en d is the time at which the shallow phase ends, for which we 
can only set a lower limit. We note that this depends weakly 
on the onset time of the shallow phase, therefore the calculated 
fiuence is correct in both presented scenarios. For comparison, 
the amount of energy released during the steep phase of the 
light curve (excluding the flare) is T ~ 2 x 10~ 8 (f sta rt/100 s)~ (X6 
erg crrT 2 . We note that the shallow phase lasts a considerable 
time. 

IZhang et all d2Q05t> propose three explanations for 
the energy injection mechanism. In the impulsive case 
dSari & Meszarosl2000l) . the central engine ejects material with 
a wide distribution of Lorentz factors. In this case, slower mov- 
ing shells will catch the fireball at a later time. We can estimate 
the minimum Lorentz factor as r m j n 2(£ , ; so 5o/no) 1 ^ 8 (l +z) 3 ^ 8 , 
where £i s0 = E^a.^ x 10 50 erg is the isotropic-equivalent en- 
ergy, and no is the external medium particle density in units 
of cm 4 . This implies that the acceleration process works 
from ultra- to mildly-re lativistic velocities. Within the putative 
Poyinting flux scenario dZhang & Kobavashil2005l) . the energy 
supply is provided by the transfer of magnetic energy to the 
fireball, and the time at which the injection stops is related 
to the ratio cr of the electromagnetic to baryonic kinetic 
energy. If this scenario is correct, we can infer a lower limit of 
o- = (fendAstart) 1 " 9 > 10-100, for q = 0.5-0, where f start < t b is 
the start time of injection. Therefore, after the end of the energy 
transfer phase, the energy of the blast-wave would be increased 
by a comparable factor. In the thi rd scenario (the p rolonged 
energy output by the central engine, Zhang & Mesza rosl200ll) . 
the end of the injection phase is simply the end of the engine 
activity. In this case, this activity produces a large amount 
of energy, particularly so since the radiative efficiency may 
be lower during the late afterglow than during the prompt 
e mission, as is generall y the case. This was previously noticed 
bv lNousek et alJ (2006)' in a sample of several Swift GRBs. 

The monitoring of XRF 050406 was discontinued 22 days 
after the trigger. By then, the source was no longer detectable 
and only a 3-cr upper limit could be drawn at » 3.6 x 10~ 4 
counts s . In order for the afterglow energy not to diverge, 
a further, late break is necessary. One interesting possibil- 
ity is that this may be due to seeing the edge of the jet. A 
steepening in the light curve is expected when the fireball 
Lorentz factor becomes comparable to the inverse of the jet 
half-opening angle. Such a late break is not unexpected for 
an XRF. The few XRFs with known redshif t dSoderberg et alJ 
l2004l iBersier et alJl2005l iFvnbo et alJl2004l) have a very low 
isotropic -energy release, and this may be at least in part ac- 
commodated if they have very wide jets. This picture is con- 
sistent with the resu lt found by Frail et al. (2001; see also 
Ghirlanda et a i1l2004 . who found that low-energy GRBs tend 
to have wider opening angles. Using the standard formalism 
dRhoadsl Fl999: Sari et aljfl999b . the jet half-opening angle is 
0j = 16 ^ 8 no /8 (^/0.2) 1/8 £ry 5 8 deg, where fj = r j;6 x 10 6 s is the 
jet break time and rj is the burst radiative efficiency. Therefore, 
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using our lower limit on the jet break time fj ^ 10 6 s, we can 
infer a lower limit on the jet half-opening angle of 16 deg. 
This value is at the high end of the distribution of jet angles 
jBloom et al.l2003bl) . 

5. Summary and conclusions 

XRF 050406 is classified as an X-ray flash, with fluence ~ 
1 x 10~ 7 erg cirT 2 (15-350 keV), a soft spectrum (T y = 2.65), 
no significant flux above ~ 50 keV and a peak energy E v < 15 
keV. Its main characteristics are however not qualitatively dif- 
ferent from those of normal GRBs. As observations accumu- 
late, it becomes clear that these two classes of phenomena share 
many properties, and both have afterglows with similar charac- 
teristics. This is a clue that both events may have a common 
origin. 

XRF 050406 is the first Swift-detected burst that showed a 
flare in its X-ray light curve, a feature now found in ~ 50 % 
of the XRT afterglows. The flare peaked at ~ 210 s after the 
BAT trigger (~ 61 s in the rest frame). The best fit of the af- 
terglow decay is obtained with a broken power law with a\ = 
1.58 + 0.17, a 2 = 0.50+j;-^, and a break at ~ 4400 s after the 
BAT trigger. The mean photon index is Fx = 2.1 + 0.3. During 
the X-ray flare a flux variation of 6F/F pea t ~ 6 in a timescale 
<Sf/fpeak <K 1 is observed, and its measured fluence in the 0.2-10 
keVbandis ~ 1.4x 10~ 8 erg cm" 2 [(2.0+ 1.4)xl0 50 erg], which 
corresponds to 1-15% of the prompt fluence. We argued that 
the flare-producing mechanism is late internal shocks, which 
implies that the central engine is still active at t ~ 210 s, though 
with a reduced power with respect to the prompt emission. We 
showed possible indications of spectral variations during the 
flare, and a flattening of the X-ray light curve after t ~ 4400 s 
in support of continued central engine activity at late times. 

Since XRF 050406 was observed, flares have been detected 
by XRT in both X-ray flashes and normal GRBs, indicating that 
flares are linked to some common properties of both kinds of 
bursts, and probably tied to their central engine. 
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Table 1. Observation log of XRF 050406. 
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Sequence 


Obs/Mode 


Start time (UT) 


End time (UT) 


Exposure" 


Time since trigger 






(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 


(s) 


(s) 


00113872000 


BAT 


2005-04-06 15:53:48 


2005-04-08 06:32:39 




-300 


00113872000 


XRT/IM 


2005-04-06 16:00:12 


2005-04-06 16:00:14 


2.5 


84 


00113872000 


XRT/PuPD 


2005-04-06 16:00:17 


2005-04-08 00:04:58 


92 


90 


00113872000 


XRT/LrPD 


2005-04-06 16:00:18 


2005-04-08 09:33:32 


681 


91 


00113872000 


XRT/WT 


2005-04-06 16:00:20 


2005-04-08 09:38:58 


3728 


92 


00113872000 


XRT/PC 


2005-04-06 16:00:26 


2005-04-08 09:49:03 


49939 


99 


00113872001 


XRT/PC 


2005-04-08 09:49:39 


2005-04-08 22:51:57 


7431 


150652 


00113872002 


XRT/PC 


2005-04-09 00:14:14 


2005-04-11 23:09:57 


12558 


202527 


00113872005 


XRT/PC 


2005-04-12 00:36:42 


2005-04-12 23:20:12 


3912 


463074 


00113872006 


XRT/PC 


2005-04-13 00:29:18 


2005-04-13 23:26:58 


11291 


549031 


00113872008 


XRT/PC 


2005-04-16 00:48:33 


2005-04-18 23:43:57 


35937 


809386 


00113872009 


XRT/PC 


2005-04-20 01:13:27 


2005-04-20 23:59:57 


17937 


1156479 


00113872010 


XRT/PC 


2005-04-21 01:19:44 


2005-04-21 23:59:58 


12546 


1243256 


00113872011 


XRT/PC 


2005-04-22 00:05:06 


2005-04-22 22:48:57 


11526 


1325178 



The exposure time is spread over several snapshots (single continuous pointing at the target) during each observation (with the exclusion 
of BAT and XRT/IM data). 



Table 2. Spectral fit results. 



Spectrum 


Photon index 


iVn 

(10 20 cm" 2 ) 


xL ( d -°- f -) 


C-stat (%) a 


Start time 
(s since T ) 


End time 
(s since T Q ) 


BAT total 

BAT T 90 

bat r 50 

BAT peak 


2 63 +042 
2 64 +0 - 46 

z " u -0.38 

2 65 +0 - 61 
2 65 +0 82 




1.4 (56) 
1.3 (56) 
1.0 (56) 
0.9 (56) 




-2.560 
-2.432 
-0.064 
0.064 


4.160 
3.648 
2.048 
1.024 


XRTWT 
XRT WT+PC 
XRTPC 
XRT PC C 


2 ir°- 31 

9 | 9+0.25 
-0.23 
9 i -2+0.44 
z " 1 -0.19 

2 06 +0 24 


2.8 b 
2.8 b 
2.8 b 
2.8 b 


1.0 (6) 
1.0 (8) 


220.7 (58.0) 
270.7 (67.1) 


92 
92 
99 
599 


596 
599 
599 
18308 



a Cash statistic (C-stat) and percentage of Monte Carlo realizations that had statistic < C-stat. We performed 10 4 simulations. 
b Fixed to the Galactic value. 

c Snapshots 1 through 4, with the exclusion of the first ~ 600 s (non piled-up data). 
Table 3. Light curve fit results' 1 . 



Model parameters Simple power law Broken power law Smoothly-joined power laws Broken power law+Gaussian 



excluding flare excluding flare excluding flare full data set 



an 1.41+ , 22 1.58*»2 1.73*$ 1.58 ±0.17 

f b (s) - (4.19^) xlO 3 (3.61+J3f) x 10 3 (4.36^ 23 )x 10 3 

a 2 0.50 ±0.14 0.42+0;!' 0.50^;| 3 

Gaussian centre (s) - - - 21 1.1 +4 4 

Gaussian width (s) - - - l 7 -9! 4 2 6 3 

X 2 mi 4.32 1.20 1.29 1.58 

d'o.f. 12 10 10 17 



We follow the notation F(t) oc r a . 



